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Kinematic Structure in the Galactic Halo at the North 
Galactic Pole: RR Lyrae and BHB Stars show different 
kinematics. 
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ABSTRACT 



Radial velocities and proper motions (derived from the GSC-II database) are 
given for 38 RR Lyrae (RRL) stars and 79 blue horizontal branch (BHB) stars in a 
~ 200 deg 2 area around the North Galactic Pole (NGP). Both heliocentric (UVW) 
and galactocentric (Vr, V^, V z ) space motions are derived for these stars using a 
homogeneous distance scale consistent with (m — M)o=18.52 for the LMC 

An analysis of the 26 RRL and 52 BHB stars whose height (Z) above the plane 
is less than 8 kpc, shows that this halo sample is not homogeneous. Our BHB sample 
(like that of Sirko et al. 2004b) has a zero galactic rotation (V^) and roughly isotropic 
velocity dispersions. The RRL sample shows a definite retrograde rotation (V^ = 
— 95±29 km s _1 ) and non-isotropic velocity dispersions. The combined BHB and RRL 
sample has a retrograde galactic rotation (V) that is similar to that found by Majewski 
(1992) for his sample of subdwarfs in SA 57. The velocity dispersion of the RRL stars 
that have a positive W motion is significantly smaller than the dispersion of those 
"streaming down" with a negative W. Also, the ratio of RRL to BHB stars is smaller 
for the sample that has positive W. 

Our halo sample occupies 10.4 cubic kpc at a mean height of 5 kpc above the 
Galactic plane. In this volume, one component (rich in RR Lyrae stars) shows retro- 
grade rotation and the streaming motion that we associate with the accretion process. 
The other component (traced by the BHB stars) shows essentially no rotation and 
less evidence of streaming. These two components have HB morphologies that suggest 
that they may be the field star equivalents of the young and old halo globular clusters 
respectively. Clearly, it is very desirable to use more than one tracer in any kinematic 
analysis of the halo. 

Key words: Stars: kinematics; Stars: horizontal branch; Stars: variables: RR Lyrae; 
Galaxy: structure; Galaxy: halo 



1 INTRODUCTION 

A well known view of the globular-cluster halo is that it con- 
sists of an Old Halo and a Young Halo whose clusters have 
different HB morphologies that can be interpreted as an age 
difference (Zinn 1993). The two systems are to some extent 
coterminous but the clusters whose R 9a ; < 6 kpc all belong 

* e-mail: kinman@noao.edu 

f NOAO is operated by AURA Inc. under contract with the Na- 
tional Science Foundation. 



to the Old Halo, while the Young Halo predominates in the 
outer parts of our Galaxy. Recent analyses of cluster ages 
(Rosenberg et al. 1999; De Angeli et al. 2005) show that the 
most metal-poor clusters are all older and coeval while the 
more metal-rich are younger and have a significant disper- 
sion in age. The outer parts of the Galaxy show a somewhat 
larger dispersion in cluster age. The Old Halo system has a 
prograde galactic rotation while that of the Young Halo is 
retrograde. The Old Halo clusters have a high ratio of blue 
horizontal branch (BHB) to RR Lyrae (RRL) stars. Layden 
(1996) found that the field RRL stars with 3 ^ K ga i ^ 6 
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kpc were not kinematically as cool as the Old Halo clusters 
since they have a larger velocity dispersion and a milder 
prograde galactic rotation. Layden also analysed the radial 
velocities of the RRL stars with |Z| > 5 kpc and whose R 9a ; 
projected on the plane were between 6 and 10 kpc; these 
showed a marginally retrograde rotation. He accounted for 
these results in terms of the preferential destruction of the 
redder-HB clusters in the inner halo and by these redder-HB 
clusters having a greater fraction of RRL stars per unit lu- 
minosity. Borkova and Marsakov (2003) analyzed the space- 
motions of local RRL stars and also found evidence for two 
metal-poor ([Fe/H]< — 1.0) systems. One (corresponding to 
the Old Halo) is associated with blue-HB clusters and has a 
spherical distribution and a slightly prograde galactic rota- 
tion. The other (corresponding to the Young Halo) is more 
spatially extended, is associated with red-HB clusters and 
has stars on eccentric and retrograde orbits. A correlation 
between the kinematics of globular clusters and their Oost- 
erhoff type was pointed out by van den Bergh (1993) and 
also by Lee & Carney (1999). They found that Oosterhoff 
type II clusters belong to an older system with prograde or- 
bits while Oosterhoff type I clusters to a younger (accreted) 
system on retrograde orbits. Lee & Carney associated the 
field RRL stars with |Z| < 3 kpc and |Z| > 5 kpc with Oost- 
erhoff types II and I respectively. The discovery of the dis- 
integrating Sgr dwarf galaxy (Ibata, Gilmore & Irwin 1994) 
gave very strong support to the idea that the outer field-star 
halo parallels the Young cluster halo and is composed of the 
tidal debris from such accreted satellites. This has prompted 
searches for other satellites and consequent non-uniformities 
in the structure of the halo. 

Surveys to discover such substructure include among 
others, a K-giant survey (Morrison et al. 2000), an APM 
carbon star survey (Totten & Irwin 1998; Totten, Irwin & 
Whitelock 2000; Ibata et al. 2001), and the Sloan Digital Sky 
Survey (SDSS) for A stars (Yanny et al. 2000) and RRL stars 
(Ivezic et al. 2000). All have shown inhomogeneities in the 
distribution of stars in the outer halo that are mostly debris 
from the Sgr dSph (Newberg et al. 2002). More recent work 
with SDSS data has shown that these "streams" and over- 
densities are very extensive (Lupton et al. 2005, Belokurov 
et al. 2006). Many of these surveys study the distant halo: 
88% of the large sample of BHB stars discovered by Sirko 
et al. (2004a) in the SDSS lie beyond 10 kpc and only 4 of 
these stars lie within 6 kpc. Thus, since proper motions are 
largely lacking for these stars, what we know of their kine- 
matics can only come from analyzing their radial velocities 
(Sirko et al. 2004b; Thorn et al. 2005). 

The kinematics of the local halo is much better known 
because accurate proper motions are available (Martin & 
Morrison 1998; Maintz & de Boer 2005) Q The mean halo 
rotation in the solar neighbourhood is close to zero (i.e. 
(V) ~ -220 km s _1 ) or is slightly prograde, and it is well 
established that the velocity dispersions in the local halo are 
non-isotropic ( o~u > <rv > o~w)- Helmi et al. (1999) studied 
97 metal-deficient ([Fe/H]^ — 1.6) red-giants and RRL stars 
within 1 kpc of the Sun and found that about ~10% of these 
had space-motions that showed that they came from a sin- 
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gle disrupted satellite. This has been confirmed with a larger 
sample (234 stars within 2.5 kpc of the Sun) by Kepley et 
al. (2006) who also found evidence for two other possible 
"streams" . 

One of the first indications of structure in the halo was 
the discovery of retrograde motion for a group of subdwarf 
halo stars in Selected Area (SA) 57 at the North Galactic 
Pole (NGP) by Majewski (1992). Later, Majewski, Munn & 
Hawley (1996) found that these halo stars were clumped in 
phase-space and [Fe/H]; in particular, they found a retro- 
grade group «V) = -275 ± 16 kins -1 ) with [Fe/H] < -0.8 
moving towards the galactic plane. Kinman et al. (1996) 
also found a preponderance of negative radial velocities (i.e. 
motion towards the plane) among BHB and RRL stars in 
the same area of the sky. A detailed analysis of the Sgr 
stream of debris using 2MASS M-giants (Majewski et al. 
2003; Law, Johnson & Majewski 2005) shows that the Sgr 
stream may be inflowing from the general direction of the 
NGP towards the solar neighbourhood. Martmez-Delgardo 
et al. (2006) have identified the Sgr leading arm with the 
overdensity in Virgo and predict highly negative radial ve- 
locities in the North Galactic Cap but not at the area of the 
Pole itself. On the other hand, neither Carney (1999) nor 
Chiba & Beers (2000) find retrograde motion among differ- 
ent samples of nearby halo stars whose Z maa; ^ 4 kpc. Val- 
lenari et al. (2006) used the Padova Galaxy model to analyse 
the proper motions in a 26.7 deg 2 near the NGP. Their halo 
component (for Z < 7 kpc) showed no significant rotation 
and had non-isotropic velocity dispersions. The colour dis- 
tributions of their halo stars suggest that they are primarily 
redder ((£> — V) ^ +0.5) giant stars. Other reports of nearby 
halo substructure include Chiba & Mizutani (2004), Alt- 
mann, Catelan & Zoccali (2005) and Meza et al. (2005) who 
discuss the debris from the globular cluster u Cen (prob- 
ably the nucleus of a disrupted dwarf galaxy). The orbit 
discussed by Chiba & Mizutani does not rise more than 4 
kpc above the plane in the direction of the NGP and so 
this debris is probably not connected with the substructure 
observed by Majewski et al. (1996) but might account for 
a halo overdensity of F and G stars observed by Gilmore, 
Wyse & Norris (2002). Duffau et al. (2006) have found a 
clump of both RRL and BHB stars at about 20 kpc that is 
coincident with an overdensity of F-type stars that were pre- 
viously found by Newberg et al. (2002). Juric et al. (2006) 
have found a significant overdensity of halo stars at 7 to 
8 kpc in Virgo although a search by Brown et al. (2004) 
for overdensities of BHB stars out to 9 kpc gave negative 
results. Clewley & Kinman (2006) have identified two new 
clumps that contain both RRL and BHB stars at distances 
of 8 and 9 kpc. There is therefore considerable evidence that 
structure exists in the nearby halo but accurate proper mo- 
tions and radial velocities are needed in order to establish 
its significance. 

This paper (a continuation of Kinman et al. 1996, 2003, 
2005) derives space motions for 117 BHB and RRL stars in 
~200 deg 2 in the direction of the NGP. The sample of BHB 
stars in SA 57 that comes from the Case Low-Dispersion Sky 
Survey (Sanduleak 1988) and that was discussed by Kinman, 
Suntzeff & Kraft (1994) has been expanded to cover about 
six Palomar Schmidt Sky Survey fields (Fig. 1) using sources 
discussed in Sec. 2.1. The RRL stars in the same area of the 
sky have also been included and in many cases re-observed. 
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A data-base containing full details of these stars is in prepa- 
ration (Kinman et al. 2007). BHB and RRL stars continue 
to be recognized as prime halo tracers; recent surveys for 
them include Sirko et al. (2004a), Vivas et al. (2001); Clew- 
ley et al. (2002, 2004, 2005); Brown et al. (2003, 2004, 2005); 
Christlieb, Beers & Thorn (2005). Red giants, subgiants and 
subdwarfs are more broadly representive of the halo but can 
be more difficult to separate from their disk counterparts 
and estimates of their distances are commonly less exact 
than those of the HB stars. A major reason for studying 
the halo in the galactic polar regions is that the confusion 
between disk and halo stars is minimized in these direc- 
tions; this is particularly important for the RRL stars which 
have a significant disk population. The galactic extinction 
is also minimal. Further, Kepley et al. (2006) have shown 
that streams may be more easily identified by their radial 
velocities at the poles. Much previous work has depended 
on analyses of radial velocities alone. Surveys for distant 
halo stars tend to be made at high galactic latitudes well 
away from the apex and antapex of galactic rotation; the 
galactic rotation component of the radial velocity is small 
compared with the random velocity components in these di- 
rections. Consequently halo galactic rotations are not well 
determined from the radial velocity alone. At the galactic 
poles, the relation between the galactic velocity vectors U, V 
and W of a star and its distance, radial velocity and proper 
motion is particularly simple. The U and V motions de- 
pend essentially only on the distance and proper motions 
while the W motion depends almost entirely on the radial 
velocity. In principle therefore, an analysis of the errors is 
more straightforward in the polar regions (and also along 
the prime galactic meridian) than in other directions. 

It follows that it is not only important to get the best 
possible distances, radial velocities and proper motions, but 
have a realistic assessment of their uncertainties. We feel 
that the proper motions that we have used (derived from 
the GSC-II database) are the best currently available. We 
have estimated their errors from the proper motions of QSO 
in the same fields. Presumably, the best proper motions will 
come from astrometric space missions such as GAIA and 
SIM (see e.g. GAIA, Concept and Technology Study Report 
2004, Table A. 4), but it will be several more years before 
these are available. The majority of our radial velocities are 
of high quality (±10 km/s); there are some of lower accu- 
racy but these do not significantly affect the conclusions of 
this study. Finally, the database itself is not complete and 
it is hoped that it can be continually updated. In particu- 
lar, it would be valuable to include the halo red giants in 
this region (Majewski 2005). Clearly this study is an ongo- 
ing project rather than one where definitive results can be 
obtained with the current data. 

We present our sample of BHB and RRL stars and the 
data used for this study in Sect. 2; we discuss our analysis 
and results in Sect. 3, and Sect. 4 contains the summary and 
our conclusions. 



2 THE DATA 

2.1 Target Selection 

Our sample of halo tracers consists of 79 BHB and 38 RRL 
stars at distances between 1.5 and 16 kpc located in an area 




Figure 1. The POSS-I fields 266, 267, 268 (top) and 320, 321, 
322 and 323 (bottom), used as first epoch plates to derive the 
proper motions of our target stars at the NGP. Right ascension, 
along the X axis, is expressed in arc degrees. The cross shows the 
position of the NGP. For comparison, the J2000 coordinates of 
the SA 57 are 13:08:38 +29:23:06 (197.158 +29.385 in decimal 
degrees). 



of approximately 22° x 12° near the NGP. This area, shown 
in Fig. [T] is the combination of 7 POSS-I fields where proper 
motions from the GSC-II database were measured. 

The BHB stars of our sample were selected from the sur- 
veys of Sanduleak (1988), Beers et al. (1996), MacConnell et 
al. (1993), and Kinman et al. (1994). They are listed in Table 
111! with a simplified nomenclature: e.g we use 16549-51 for 
BPS BS 16549 0051 (Beers et al. 1996), AF-003 for Case A-F 
003 (Sanduleak 1988), and SA57-001 for SA 57 001 (Kinman 
et al. 1994). The single ID used in Table [PT1 is that of the 
earliest of these catalogues in which the star is mentioned; 
alternative ID will be given in Kinman et al. (2007). These 
candidate BHB stars were confirmed by uBV-photometry 
(Kinman et al. 1994) and spectroscopy. Most of the con- 
firming spectra were taken at the Kitt Peak 4-m Mayall 
telescope (Kinman et al. 1996). 

The RRL stars (Table I12|l were mostly taken from 
the General Catalogue of Variable Stars (GCVS, Kholopov 
1985) and subsequent Name-Lists, and have the traditional 
identification by constellation; those without GCVS names 
are taken from Kinman (2002). These are NSV5476, AF- 
791, AF-155, SA57-19, SA57-47, SA57-60 and AF-882. The 
data for AF-031 and AF-042 are as yet unpublished. 

Intensity-weighted mean V magnitudes of the RRL 
stars are derived from recently observed light curves. Spec- 
tra of several RRL stars were obtained at the 3.5-m TNG 
telescope. Only 6 of these RRL stars and none of the BHB 
stars in our current sample are included in the catalogue of 
halo stars by Beers et al. (2000). 

K-magnitudes are available in the 2MASS catalogue for 
all our BHB and RRL stars. The data from which the kine- 
matic properties of these stars were derived are given for the 
BHB and RRL stars in Tables [TT1 and [121 respectively. More 
detailed descriptions of this photometric and spectroscopic 
data will be given in Kinman et al. (2007). 
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2.2 Reddenings, Absolute Magnitudes and 
Distances 

2.2.1 Reddenings 

The reddening E(B — V) was taken from Schlegel et al. 
(1998). The reddening corrections for other colours are 
derived from the relations E(V - K )=2.75E(B - V), 
A V =3.1E(B - V) and A K =0.35E(B - V) (Cardelli et al. 
1989). 

2.2.2 Absolute Magnitudes and Distances: RRL stars. 

We consider three ways to estimate the absolute magnitudes 
and hence parallaxes of our RRL stars: 

(i) The visual absolute magnitude My is given in terms 
of the metallicity [Fe/H] by the empirical relation: 

M v = Q.2U[Fe/H\ + 0.86 (1) 

We adopt the coefficients given by Clementini et al. (2003) 
which imply that M v = 0.54 for [Fe/H] = -1.5. They give 
an LMC modulus of ~18.52 using the (Vo) = 19.064 that 
Clementini et al. (2003) derive from their sample of about 
100 RRL stars in the LMC bar. This Mv estimate is there- 
fore squarely in the middle of the range derived from recent 
LMC distance determinations, most notably 18.50 (Freed- 
man et al. 2001) and 18.54 (Tammann et al. 2003). The 
metallicity, [Fe/H], however, is known for only just over half 
of our RRL sample (see Table [T2")). We assume [Fe/H] = -1.6 
for the remaining stars since this is close to the mean metal- 
licity of a halo population. An error of ±0.5 dex in [Fe/H] 
leads to an error of ±0.1 mag. in the distance modulus and 
about 5% in the parallax. The parallaxes derived by this 
method (HV) and their errors (an) are given in columns 8 
and 9 respectively in Table IT21 

(ii) The infrared absolute magnitude Mk is derived from 
the metallicity and period. In the case of the type ab RRL 
stars we use the relation in the form given by Nemec et al. 
(1994): 

M K = -2.40 log P + 0.06\Fe/H] - 1.06 (2) 

where the zero point is on the same scale as that given in 
(iii) below. In the case of the c-type variables (indicated by 
an asterisk against their logP in column 4 of Table [T2")) . 
the periods must be "fundamentalized" before using them 
in this relation. For this we have assumed that the ratio of 
the first overtone (c-type) to fundamental (ab-type) period 
is 0.745 (Clement et al. 2001). The parallaxes derived in this 
way (II k ) are given in column 11 of Table [121 

(iii) The infrared absolute magnitude Mk is derived from 
the (V — K)o colour using an empirical relation derived from 
data for non-variable stars in the globular clusters M3 and 
M13 (Valenti et al. 2004). These two clusters have similar 
metallicities which are both close to the mean for the field 
halo population. We assumed reddenings E(£> — V) of 0.01 
and 0.02 and distance moduli (m— M)o of 15.07 and 14.25 for 
M3 and M13 respectively. These M3 and M13 data (Figure 
[2| can be fitted by the 3-a rejection cubic polynomial: 

M K = 0.878 - 1.812(V - K) ± 0.675(V - K)l 

-0.183(V - K)% (3) 

This relation only holds for stars with the metallicity of 



M3 and M13. The M K vs. (V - K) Q relations for ZAHB 
stars of different [Fe/H] were obtained from the models of 
VandenBergh et al. (2000). These models were kindly trans- 
formed into Mk vs. (V — K)o by A. Sollima using the 
color-temperature relation given by Bessell et al. (1998). 
The ZAHB relations show that in the colour range of the 
RRL stars an expression for Mk needs a metallicity term ~ 
+0.18[Fe/H]. Equation (3) is therefore amended to: 

M K = 1.166 ± 0.l8[Fe/H] - 1.812(V - K)o 

±0.675(V - Kf - 0.183(V - K)% (4) 

The parallaxes (JIhb), derived from equation (4), are 
given in column 12 of Table [121 

The mean difference between the Mk derived from 
equation (4) and that from equation (2) for our 38 RRL stars 
is ±0.025±0.022 mag with a dispersion of 0.134 mag; this is 
satisfactorily small. Mean infrared magnitudes < K > were 
derived from the single epoch 2MASS magnitudes using the 
template K light curves given by Jones et al. (1996) and the 
most appropriate ephemerides available. The if-amplitudes 
are such that errors of ~0.20 mag. in < K > are possible if 
the estimated phases are in error by 0.1. The letter following 
< K > in Table [T2] gives the photometric quality Q of the 
2MASS data. Thirty five of the RRL sample have quality A 
(o-k <0.10 mag); for these TLk /Tlv = 1.010 ± 0.010 with 
a dispersion of 0.062. The twenty with no metallicity esti- 
mate (for which we assumed [Fe/H] — —1.6) have Hk/Rv 
= 1.022± 0.014 and a dispersion of 0.064 while for the fif- 
teen variables for which [Fe/H] is known, this ratio is 0.996± 
0.015 with a dispersion of 0.058. The parallaxes derived from 
infrared magnitudes provide a useful confirmation of the va- 
lidity of the parallax (Hv) derived from visual magnitudes, 
but the greater uncertainty in the infrared magnitudes com- 
pared with the visual magnitudes of our program stars and 
the greater uncertainties of equations (2), (3) and (4) com- 
pared with equation (1) has led us to adopt (ILv) alone for 
the calculation of the RRL distances. 

2.2.3 Absolute Magnitudes and Distances: BHB stars. 

For the BHB stars, My can be estimated from the Mv vs. 
(B — V) relation given by Preston et al. (1991). We adjusted 
this relation so as to be consistent with the absolute scale 
that we have adopted for the RRL stars, to get: 

M v = 1.00 - 4.423(5 - V)o + 17.74(5 - V)l 

-35.73(5 - V)o (5) 

The parallax derived from equation (5) is given as 
n_f/ sp in Table 1111 Various problems involved with deriv- 
ing Mv from (B — V)o have been discussed by Brown et 
al. (2005). One concern is that the data for the 15 globular 
clusters on which this relation is based are now quite old. 
For comparison, we therefore derived a similar cubic from 
more recent data for the intermediate-metallicity clusters 
M3 (Ferraro et al. 1997) and M13 (Paltrinieri et al. 1998): 

2 The letters A, B, C, D and U correspond to the observation 
having a S/N of >10, >7, >5, <5 and being an upper limit re- 
spectively. 
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Table 1. Mean galactocentric radial velocities (RV 9a ;) for our 
sample of NGP halo stars at different heights (Z) above the plane 
and Right Ascensions (RA). 
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Figure 2. The combined HB of the globular clusters M3 and M13 
in the IR plane Mjf vs (V-K)o (data from Valenti et al. 2004). 
The data have been corrected for reddening and shifted according 
to the individual distance moduli, as specified. The line shows the 
relation in Eq. (3). 



M v = 1.145 - 5.615(5 



V)o + 16.28(5 
-14.93(B 



V)l 
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The parallax derived from equation (6) is given 
as Hhbc in Table 1111 The mean value of the ratio 
U-hbc/IIhbp is 1.032±0.003 with a dispersion of 0.027. This 
ratio is adjusted to be unity at the blue edge of the instabil- 
ity strip and increases to 1.107 for the bluest BHB star in 
our sample ((£> — V)o = —0.08). We assume that the differ- 
ence (A) between these parallaxes is a measure of the likely 
systematic error in IIhbp- If o is the error in IIhbp (as 
deduced from the error in (B — V)o) then its adopted error 
(an given in column 9 of Table [11} is \/(A 2 + a 2 ). 

Equation (3) can also be used to derive a parallax for 
the BHB stars (JIhvk in column 10 of Table [TTJ . Fifty of 
them have 2MASS K magnitudes of quality A; for these 
the ratio Hhvk /IIh bp = 0.950±0.005 with a dispersion of 
0.036. This ratio (adjusted to be unity at the blue edge of 
the instability strip) decreases to 0.888 for the bluest star in 
this sample ((B — V)o = —0.02). Equation (3) is least-defined 
for its bluest colours since these correspond to the faintest 
HB stars in the calibrating globular clusters. The derivation 
of parallaxes from (V — K)o for these stars is likely to have 
significant advantages over the use of (B — V)q, but currently 
we need better calibration and more accurate magnitudes 
for the fainter program stars before the method can be used 
with confidence. Consequently only Hhbp is used in this 
paper. 
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Figure 3. The distributions of galactocentric radial velocity as a 
function of Z distance above the galactic plane, and at different 
Right Ascensions (RA). 



2.3 Radial velocities 

The radial velocities (RV) for our halo stars were either ob- 
tained from spectra obtained at the 4m-RC (KPNO) and 
3.5m-LRS (TNG) spectrographs or from the literature. A 
more detailed discussion that includes the sources of these 
velocities and the individual galactocentric velocities will be 
given in Kinman et al. (2007) . The heliocentric RV (with er- 
rors) of the program stars are listed in Tables [TT1 and [121 We 
have no RV for 7 of the RRL stars; these all have Z > 9 kpc. 



In this analysis, we determine galactocentric ra- 
dial velocities (HV ga l), by assuming a solar motion of 
(U,V,W) Q = (10.0, 5.25, 7.17) kms" 1 with respect to the 
LSR (Dehnen & Binney 1998) and an LSR velocity of 220 
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km s" 1 . The mean galactocentric radial velocities and their 
corrected dispersion^ are given in Table Q] their distribu- 
tions are shown in Fig. [3] We divide our sample of stars into 
those that lie below and above Z = 4 kpc; the group with Z 
< 4 kpc does have a smaller dispersion, but the difference is 
only significant at the 95% level. There are eight stars in the 
Z < 4 kpc sample that have a (RVgai) of +14.8 kms -1 and 
a dispersion of only 4.2 kms -1 ; a W-test (Shapiro & Wilk 
1965) of the whole Z<4 kpc sample, however, shows no de- 
parture from a Gaussian distribution. The dispersions found 
for the sample with Z > 4 kpc are in general agreement with 
those found for halo stars by Clewley et al. (2004) i.e. 108 
± 10 kms -1 at distances 11-52 kpc from the Sun, and by 
Sirko et al. (2004a) i.e. 99.4±4.3 kms -1 for their stars out 
to 30 kpc. They are, however, smaller than the dispersion of 
120 kms" 1 (out to 30 kpc) found by Battaglia et al. (2005). 

The velocity distribution should be Gaussian with a 
zero mean velocity (Harding et al. 2001) but our Z > 4 kpc 
sample has a mean galactocentric RV of — 26.1±10.7 kms" 1 
and an excess of large negative velocities. This excess is most 
pronounced for the 23 stars whose R.A. > 13 h and whose 
mean galactocentric RV is — 49±21 kms -1 . Individually the 
two samples with RA greater and less than 13 :00 m show 
no significant departure from a Gaussian distribution on a 
W-test (Shapiro & Wilk 1965). 

2.4 Proper Motions and their Errors 

The proper motions are derived from the plate material used 
for the construction of the GSC-II catalogue (Lasker et al. 
1995; McLean et al. 2000). We used multi-epoch positions 
that were derived from digitized Schmidt plates from the 
POSS-I, Quick-V, and POSS-II surveys (see Tables and 
[3]) . These cover a time-baseline of about 40 years and should 
allow us to get a precision of a few mas per year. We used 
the procedure described by Spagna et al. (1996) to get the 
relative proper motions. These were then transformed to the 
absolute reference frame by forcing the extended extragalac- 
tic sources in the field to have no tangential motion. The 
expected zero point accuracy is ^ 1 mas yr -1 . The proper 
motions and their individual errors for our program objects 
are listed in Tables [TT1 and [121 their formal r.m.s. errors are 
^ 3 masyr -1 . Although these proper motions derived from 
the GSC-II are among the most accurate that are available 
for the large-field surveys, since an error of 3 mas yr -1 corre- 
sponds to 14.2 km s -1 at 1 kpc, the errors in the computed 
space motions of halo stars from their proper motions alone 
will be comparable with the space motions themselves at a 
distance of ~10 kpc. 



Proper motions for six of our brightest program stars 
are also given in the Hipparcos and Tycho catalogues (avail- 
able in the NOMAD catalogue, Zacharias et al. 2004). These 
proper motions are compared in Table U where it is seen 
that the differences are largely within the expected errors. 
These bright stars have large images on our survey plates 

4 The method of correction is given in Sec. 3.1. 



Table 2. Field centers of the POSS-I and Quick-V Schmidt plates 
(6.4° X 6.4°) used in the proper motion determinations. 



Field 


RA 


DEC 


Epoch 


Epoch 


Epoch 




(deg) 


(deg) 


POSS-I O a 


POSS-I E b 


Quick-V c 


266 


183.82 


35.19 


1956.350 


1956.349 


1983.133 


267 


190.76 


35.20 


1950.365 


1950.364 


1983.133 


268 


197.69 


35.22 


1950.367 


1950.368 


1983.133 


320 


183.83 


29.19 


1955.293 


1955.293 


1983.293 


321 


190.28 


29.20 


1950.275 


1950.275 


1983.293 


322 


196.73 


29.22 


1955.288 


1955.287 


1982.386 


323 


203.18 


29.25 


1950.428 


1950.425 


1982.387 



a Photographic O band (103aO unfiltered). 

6 Photographic E band (103aE+red Plexiglass 2444). 

c Photographic Via band (IlaD + Wratten 12). 



Table 3. Field centers of the POSS-II Schmidt plates (6.4° x6.4°) 
that were used for the second epoch positions in the proper motion 
determinations. 



Field 


RA 


DEC 


Epoch 


Epoch 


Epoch 




(deg) 


(deg) 


POSS-II J a 


POSS-II F b 


POSS-II N c 


379 


180.64 


34.72 


1996.216 


1989.042 


1998.373 


380 


186.62 


34.72 


1991.184 


1991.184 


1995.381 


381 


192.60 


34.73 


1990.211 


1990.211 


1995.384 


382 


198.58 


34.74 


1988.279 


1989.026 


1990.343 


441 


184.63 


29.72 


1993.285 


1993.217 


1996.312 


442 


190.36 


29.73 


1990.225 


1990.225 


1990.080 


443 


196.10 


29.73 


1995.233 


1995.233 


1993.288 


444 


201.83 


29.74 


1995.307 


1995.307 


1998.376 



a Photographic Bj band (IIIaJ+GG385). 
b Photographic R F band (IIIaF+RG610). 
c Photographic I N band (IV-N +RG9). 



and probably have the poorest proper motions in our survey. 
Consequently this comparison will give limited information 
on the systematic errors that may be present in our whole 
survey. It is therefore crucial to have another independent 
estimate of the size of the proper motion errors. We there- 
fore selected some 74 QSO brighter than 18th magnitude 
that were listed in our fields by Hewitt & Burbidge (1993) 
and Veron-Cetty & Veron (2001) and 62 type C (compact) 
objects from the Kiso Schmidt Survey for UV-excess galaxies 
(Takase & Miyauchi-Isobe 1993). These compact UV-excess 
galaxies have sufficiently "stellar" images for them to be 
suitable for use as positional standards. As we noted pre- 
viously (Kinman et al. 2003), some of these QSO showed 
suprisingly large proper motions (^10 mas yr -1 ) and spec- 
tra (kindly taken by Arjun Dey and Buell Jannuzi with the 
Kitt Peak 4.0-m telescope) showed that five of them were 
stars. We therefore rejected 15 of these objects whose proper 
motions (in either coordinate) were more than 10 mas yr -1 . 
Table [5] gives the mean proper motion in each field for the 
remaining 121 objects. The fields are defined by the ID num- 
ber of the POSS-II and POSS-I plates and some fields are 
labelled with an "a" or "b" in order to identify partially 
overlapping POSS plates. Objects that fall on overlapping 
fields have multiple measures; these have been considered 
as separate measures. The unweighted mean of these 122 



Kinematic Structure in the Galactic Halo at the North Galactic Pole. 7 



Table 4. Comparison of our proper motions (O) with those from the Hipparcos (H) and Tycho (T) proper motions as given in the 
NOMAD catalogue Zacharias et al. (2004). 



ID 


Ha 


m 


Source 


Ma 


W 


Source 


Aflat 


A /is 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


AF-113 


-25.9±1.9 


-35.8+1.7 


O 


-26.3±0.9 


-33.8±0.8 


T 


+0.4±2.1 


-2.0±1.9 


SA57-007 


-13.9±1.6 


-28.5±2.2 





-13.1±0.8 


-24.1±1.1 


T 


-0.8±1.8 


-4.4±2.5 


AF-848 


-09.8±3.0 


-18.0±2.8 





-11.8+0.7 


-20.5±0.9 


T 


+2.0±3.1 


+2.5±2.9 


AF-914 


-36.2+3.4 


-15.0±3.4 





-30.5+0.7 


-13.1±0.7 


T 


-5.7±3.5 


-1.9±3.5 


S-Com 


-19.6±4.4 


-18.4±3.9 





-16.4±3.4 


-16.5±2.2 


H 


-3.2±5.6 


-1.9±4.5 


U-Com 


-49.7±2.7 


-16.4±2.1 





-44.4±3.6 


-13.8+2.0 


H 


-5.3+4.5 


-2.6±2.9 



Table 5. Proper motion systematic errors from QSO and Galax- 
ies. The field identification in column 1 shows the ID numbers of 
the POSS-II and POSS-I plates. Column 4 lists the number of 
QSO+Galaxies (Q+G) measures, and column 5 lists the number 
of BHB+RRL stars present in the corresponding fields. 



Field 


fi a ± err. 


fig i err - 


No. of 


No. of 




mas yr~ 


mas yr 


Q+G 


stars 


379-266 








8 


380a-267 


-3.0 ±. •• 


-2.0 + 


1 


5 


380b-266 


-2.0 ±- • ■ 


+4.0 ± • • ■ 


1 


4 


381a-267 


-1.68 ± 0.73 


-0.81 ± 1.35 


9 


11 


381b-268 


+0.27 ± 1.36 


-0.29 ± 1.15 


7 


1 


382-268 


-0.37 ± 0.87 


+0.49 ± 0.73 


23 


9 


441-320 


-0.24 ± 0.70 


+0.84 ± 0.95 


14 


28 


442-321 


+0.27 ± 0.50 


-0.66 ± 0.59 


30 


26 


443-322 


+0.16 ± 0.62 


-0.37 ± 0.69 


29 


21 


444a-322 








1 


444b-323 


-1.56 ± 0.74 


+0.78 ± 1.26 


8 


3 


ALL 


-0.25 ± 0.28 


-0.07 ± 0.31 


122 


117 


FIVE FIELDS 


-0.14 ± 0.31 


-0.14 ± 0.34 


105 


95 



independent measures is given as "ALL" in Table [5] and is 
zero within its errors. The r.m.s. uncertainty of this result 
is ~0.3 mas yr -1 in each coordinate. Some 86% of the QSO 
and compact galaxies and 81% of the program stars are in 
the five fields 381a-267, 382-268, 441-320, 442-321 and 443- 
322. The unweighted means of the proper motions for these 
five fields are also zero within their errors. While we can- 
not be certain that there are not larger systematic errors in 
individual fields, the available data indicate that the r.m.s. 
systematic error in the proper motion is ~0.3 mas yr^ 1 in 
each coordinate for the data as a whole. The extragalactic 
standards are somewhat fainter on average than the pro- 
gram stars so that we must also be aware that magnitude- 
dependent errors may also be present. 

3 ANALYSIS AND RESULTS 

3.1 The space-velocities U,V & W 

The heliocentric space- velocity components U, V, and W 
were derived from the data listed in Tables [11] and [12] We 
used the program by Johnson & Soderblom (1987) (updated 
for the J2000 reference frame and further updated with the 
transformation matrix derived from the Vol. 1 of the Hip- 



parcos data catalogue). This program gives a right-handed 
system for U, V and W in which these vectors are positive 
towards the directions of the Galactic centre, Galactic rota- 
tion and the NGP respectiveljjf]. As noted in Sect. 2.2.2 and 
2.2.3, Hhbp and HV were adopted for the parallaxes of the 
BHB and RRL stars respectively. No radial velocities are yet 
available for seven of the RRL stars. We therefore assumed, 
when calculating their U,V & W velocities, that their radial 
velocity is zero with an error of 150 kms -1 . This should give 
acceptable U and V velocities since these depend almost en- 
tirely on the proper motions; their W velocities must, of 
course, be discarded. These values of U,V & W together 
with the height (Z) of the star above the Galactic plane and 
the galactocentric distance (R ga i) are given in Table [6] and 
Table [7] for the BHB and RRL stars respectively. We assume 
a solar galactocentric distance of 8.0 kpc. 



Table [8] gives the mean values <U>, <V> and <W> 
of these space- velocities for various ranges of Z. The mean 
values and dispersion given in parentheses were obtained by 
trimming the 10% of the sample that have the most extreme 
values. This has little effect on the mean values and will not 
be discussed further. The rms dispersions in these quantities 
must be corrected for the errors in the individual U,V & W 
space- velocities (given in Tables [6] and [7} . Following Jones 
& Walker (1988), if the observed dispersion in U is Disp(U), 
and £ i is the error in U of star i, then the corrected dispersion 
a u is given by: 

n 

al = {Disp{U)f--Y^ei (7) 

i=l 

and similarly for the space-velocities V and W. These cor- 
rected dispersions are given in columns (5), (7) and (10) 
of Table [5] For a given proper motion error, the space- 
velocity error increases with distance; when it becomes 
comparable with the intrinsic velocity dispersion (at Z~10 
kpc), the correction becomes unreliable. Table [5] also gives 
<U>,<V> & <W> and their dispersions for three halo 
samples from the solar neighbourhood. Local A is a sample 
of stars that were confirmed as BHB stars by high-resolution 
spectroscopy (Kinman et al. 2000) ; these stars have Hippar- 
cos proper motions and radial velocities with errors of a 

5 Other authors (e.g. Majewski 1992) use a left-handed system 
in which U is positive towards the Galactic anticentre. 
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Table 6. The BHB stars of our sample. The U,V and W velocities (in km s 1 ) are heliocentric. The velocities V r , and N z (in km 
s _1 ) arc in galactoccntric cylindrical coordinates. For further details, see text. 



ID 


Z (kpc) 


R„ ai (kpc) 


U 


V 


W 




V 4> 




(!) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


16549-51 


3.6 


9.5 


-192±033 


-246±026 


+007±049 


-183±033 


-013±026 


-014±049 


AF-003 


14.8 


18.5 


+199±296 


-209±323 


-039±071 


+207±296 


+028±323 


+032±071 


AF-006 


3.2 


9.1 


+011±048 


-338±042 


+013±013 


+022±048 


-106±042 


-020±013 


AF-011 


10.0 


14.1 


+418±198 


-534±277 


-021±037 


+432±198 


-295±277 


+014±037 


AF-013 


15.6 


18.9 


-130±153 


-405±127 


+084±048 


-106±153 


-183±127 


-091±048 


AF-727 


4.3 


9.9 


+048±033 


-024±036 


-012±012 


+062±033 


+206±036 


+005±012 


AF-022 


12.6 


16.2 


+072±114 


-370±082 


-037±022 


+085±114 


-136±082 


+030±022 


AF-729 


8.0 


12.5 


+021±072 


-331±074 


-016±041 


+027±072 


-100±074 


+009±041 


AF-029 


6.6 


11.4 


+079±081 


-105±064 


-045±017 


+090±081 


+125±064 


+038±017 


AF-030 


5.8 


10.6 


-077±063 


-324±059 


+016±014 


-065±063 


-094±059 


-023±014 


16022-26 


5.8 


10.5 


+013±050 


-288±065 


-047±049 


+025±050 


-055±065 


+040±049 


AF-038 


8.7 


12.8 


-078±089 


-423±127 


-090±017 


-070±089 


-191±127 


+083±017 


AF-039 


12.2 


15.7 


-280±091 


-217±094 


+011±041 


-271±091 


+006±094 


-018±041 


AF-041 


7.4 


11.8 


-148±058 


-433±066 


+077±013 


-141±058 


-200±066 


-084±013 


AF-045 


3.4 


9.1 


-055±019 


-067±020 


+007±010 


-049±019 


+164±020 


-014±010 


AF-046 


11.0 


14.6 


-136±112 


-192±162 


+040±017 


-128±112 


+036±162 


-047±017 


AF-048 


3.6 


9.3 


-038±034 


-211±037 


+018±011 


-029±034 


+021±037 


-025±011 


AF-050 


15.2 


18.3 


+042±216 


-384±180 


+019±050 


+059±216 


-149±180 


-026±050 


AF-052 


6.0 


10.7 


-116±110 


-293±090 


-088±018 


-106±110 


-062±090 


+081±018 


AF-053 


6.8 


11.4 


+148±059 


-402±091 


+047±015 


+153±059 


-173±091 


-054±015 


AF-063 


12.4 


15.6 


-079±096 


-316±075 


-117±015 


-069±096 


-085±075 


+110±015 


16026-28 


3.2 


8.9 


+033±024 


-055±022 


-085±010 


+038±024 


+178±022 


+078±010 


AF-754 


7.8 


12.2 


+012±077 


-252±120 


+008±045 


+019±077 


-021±120 


-015±045 


AF-755 


7.8 


12.2 


-031±098 


-296±113 


-101±044 


-027±098 


-062±113 


+094±044 


AF-068 


5.8 


10.4 


-001±065 


-195±052 


-136±012 


+008±065 


+037±052 


+129±012 


AF-070 


6.9 


11.1 


-094±052 


-231±066 


+166±011 


-085±052 


-001±066 


-173±011 


AF-073 


16.8 


19.6 


+233±125 


-238±122 


-001±042 


+241±125 


-023±122 


-006±042 


AF-075 


8.3 


12.1 


+231±103 


-365±089 


-020±014 


+238±103 


-137±089 


+013±014 


AF-076 


9.9 


13.3 


+065±102 


-285±066 


+046±013 


+074±102 


-054±066 


-053±013 


AF-077 


9.2 


12.9 


+175±146 


-329±166 


-126±022 


+180±146 


-104±166 


+119±022 


AF-078 


4.4 


9.5 


+118±027 


-087±037 


-165±010 


+128±027 


+144±037 


+158±010 


AF-769 


3.5 


9.1 


-116±036 


-393±046 


+021±011 


-112±036 


-158±046 


-028±011 


16026-67 


2.9 


8.6 


-041±037 


-338±055 


-040±010 


-031±037 


-106±055 


+033±010 


AF-100 


3.9 


9.1 


-181±032 


-305±040 


+049±010 


-172±032 


-072±040 


-056±010 


16466-08 


6.3 


10.7 


-123±076 


-257±062 


+025±040 


-115±076 


-018±062 


-032±040 


AF-108 


3.8 


9.0 


-072±032 


-492±046 


-069±010 


-064±032 


-260±046 


+062±010 


AF-111 


5.0 


9.8 


-282±047 


-087±045 


-037±011 


-266±047 


+157±045 


+030±011 


AF-112 


2.6 


8.6 


+069±023 


-089±022 


+009±010 


+080±023 


+142±022 


-016±010 


AF-113 


1.6 


8.3 


-025±014 


-341±019 


+020±010 


-017±014 


-109±019 


-027±010 


AF-115 


8.1 


11.9 


-017±060 


-397±056 


+041±011 


-020±060 


-164±056 


-048±011 


16466-15 


2.9 


8.7 


-052±042 


-079±052 


+146±011 


-040±042 


+154±052 


-153±011 


AF-130 


11.4 


14.2 


+063±090 


-209±068 


+113±010 


+073±090 


+021±068 


-120±010 


AF-131 


9.2 


12.6 


-012±080 


-314±078 


+059±011 


-007±080 


-082±078 


-066±011 


AF-134 


2.6 


8.6 


+220±023 


-204±019 


+002±109 


+230±023 


+023±019 


-009±010 


16031-44 


4.2 


9.1 


+161±040 


-203±056 


+011±010 


+170±040 


+030±056 


-018±010 


15622-48 


5.7 


9.9 


+078±049 


-279±064 


-069±010 


+087±049 


-047±064 


+062±010 


15622-07 


2.7 


8.5 


+033±033 


-172±039 


-037±010 


+042±033 


+060±039 


+030±010 


AF-138 


5.2 


9.8 


-053±065 


-325±059 


+124±011 


-049±065 


-091±059 


-131±011 


15622-09 


3.0 


8.5 


+086±025 


-178±032 


-122±010 


+095±025 


+054±032 


+115±010 


AF-797 


8.1 


12.0 


-144±072 


-350±050 


+119±039 


-149±072 


-100±050 


-126±039 


AF-804 


4.5 


9.5 


-174±043 


-243±034 


+114±010 


-165±043 


+002±034 


-121±010 


SA57-001 


5.7 


9.9 


+228±069 


-445±060 


-113±010 


+234±069 


-216±060 


+106±010 


SA57-006 


4.4 


9.2 


-191±053 


-257±050 


+119±010 


-183±053 


-021±050 


-126±010 


SA57-007 


1.6 


8.2 


+030±013 


-240±026 


-084±049 


+039±013 


-008±026 


+077±049 


SA57-009 


11.8 


14.3 


-169±136 


-716±117 


+067±010 


-181±136 


-477±117 


-074±010 


SA57-017 


3.3 


8.6 


-239±069 


-356±063 


+004±010 


-231±069 


-122±063 


-OlliOlO 


SA57-021 


10.3 


13.2 


-399±134 


+066±146 


-085±042 


-358±134 


+336±146 


+078±042 


SA57-029 


12.6 


15.0 


-009±099 


-307±105 


-161±040 


-008±099 


-075±105 


+154±040 


SA57-032 


7.8 


11.3 


-074±069 


-406±084 


+017±012 


-078±069 


-169±084 


-024±012 


SA57-036 


11.9 


14.3 


+003±123 


-398±141 


+012±011 


+001±123 


-166±141 


-019±011 
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Table 6. continued 





1D 


Z (kpc) 


R ga i (kpc) 


U 


V 


W 


Vr 


V* 


v z 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


SA57- 


-040 


7.4 


10.8 


no o 1 n Ar\ 

— 033±049 


100 I n a c 

— 128±046 


i nor 1 I n 1 n 

+036±010 


ni n I n a n 

— 019±049 


i i n cr I n ,1 f 

+ 105±046 


n a o I m n 

— 043±010 


SA57- 


-041 


13.0 


15.1 


— 124±115 


—341+114 


i -\ a n I mi 

+146±011 


— 124+115 


nnn I 1 1 a 

— 099±114 


— 153±011 


SA57- 


■045 


8.2 


11.5 


i on/" 1 i n rr n 

+206±050 


O A o | r\ a r- 

— 242±U45 


1101 m n 

— 113±010 


i oi o I n crn 

+213±050 


non I c\ a 

— 030±045 


i 1 nc i n i n 

+ 106±010 


Al- 


-825 


4.8 


9.5 


n a n I nno 


1 n 1 I non 

— iyi±030 


ii/i I n a n 

— 114±U49 


no t i noo 

— 027±023 


i r\ A *> I non 

+043±030 


i 1 n t i n ,i n 

+ 107±049 


SA57- 


-046 


16.1 


17.7 


— Illil43 


— 483±157 


i i t a I n ,1 n 

+ 174±040 


— 123±143 


— 241±157 


1 o 1 \ r\ a r\ 

— 181±040 


SA57- 


-055 


11.6 


14.0 


i 1 o f i i/ii 
+ 136±141 


Q 1 ET 1 1 no 

— 315±102 


— 201±014 


I 1 Q A I 1 A 1 

+ 134±141 


i nn i i no 

— 100±103 


+ iy4±014 


AF- 


-841 


4.5 


9.2 


— 012±029 


— 273±035 


+006±010 


— 006±029 


— 041±035 


— 013±010 


AF- 


-848 


1.4 


8.1 


+014±019 


-144±021 


-099±010 


+025±019 


+088+021 


+092±010 


SA57- 


-066 


6.1 


10.0 


-009+087 


-168±075 


-225±013 


+005±087 


+064+075 


+218±013 


AF- 


-854 


6.9 


10.9 


-295±058 


-211±052 


+019±049 


-279±058 


+067±052 


-026±049 


SA57- 


-080 


4.6 


9.1 


-083±038 


-192±048 


-024+010 


-072+038 


+044±048 


+017±010 


AF- 


-866 


4.1 


9.0 


+030±042 


-245±050 


-003±050 


+038±042 


-016±050 


-004±050 


SA57- 


-087 


10.9 


13.4 


-077±135 


-308±135 


-074+018 


-079±135 


-064±135 


+067±018 


AF- 


-900 


3.4 


8.6 


-018±020 


-215±025 


+128±049 


-008±020 


+018±025 


-135±049 


SA57- 


-111 


13.6 


15.4 


-160±089 


-448±094 


-080±014 


-191±089 


-181±094 


+073±014 


AF- 


-909 


7.0 


10.4 


+060±042 


-214±043 


-121±011 


+071±042 


+010±043 


+114±011 


AF- 


-914 


1.7 


8.2 


-162+031 


-249±037 


+151±049 


-154±031 


-010±037 


-158±049 


AF- 


-916 


4.7 


9.2 


+073±033 


-130±036 


-115±049 


+091±033 


+094±036 


+108±049 


AF- 


-918 


4.1 


8.7 


-055±041 


-188±039 


+001±049 


-043±041 


+047±039 


-008±049 



few km s ; their parallaxes and space- velocities were de- 
termined in exactly the same way as for our program BHB 
stars. Local B is the HALOl sample of local RRL stars (with 
[Fe/H] < -1.3) from Martin & Morrison (1998). Local C 
is their HAL02 sample in which the total space-velocity 
(which is defined by equation (8) in Sec. 3.4) was used to 
exclude thick disk stars (as shown in our Fig. |6}. Martin 
& Morrison used M v = +0.73 (at [Fe/H] = - 1.9) to de- 
termine the distances of their RRL stars; they plot the U, 
V & W and their dispersions for their HAL02 sample as a 
function of My in Figs. 6 & 7 respectively of their paper. 
We used these plots to derive the space-motions of their 
HAL02 sample for our assumed My of +0.54 at [Fe/H] = 
— 1.5; these are given as Local D. Changes in Mv have the 
greatest effect on the V space-velocity and the U velocity 
dispersion in these local halo samples. 

If we assume that <Vlsr > = —220 km s _1 for zero 
halo rotation, the corresponding < Vhei > will be ~— 225 
km s^ 1 (see Sect. 2.3). Older local halo samples (see Table 
5 of Martin & Morrison) show slightly prograde <V>. It 
is not clear how much these prograde <V> are a result of 
the adopted distance scales and/or the presence of disk stars 
in these halo samples. Most local samples have the problem 
that either (a) they may contain some thick-disk stars (as 
in the case of the HALOl) or (b) may have some kinematic 
bias which is introduced in the attempt to remove the disk 
component (e.g. HAL02) or which is present in their original 
selection (e.g. Local A). The differences between the space- 
motions of the different local samples suggests that these 
putative systematic sampling errors are comparable in size 
to their statistical errors. 

3.2 The Halo Rotation <V> 

We see in Table[8]that the mean value <V> of our halo sam- 
ple becomes increasingly negative (corresponding to more 
retrograde orbits) with increasing height Z above the galac- 
tic plane. This change in <V> is several times its calculated 



error and so is formally significant. The V of the individual 
stars are shown as a function of Z in Fig. |4jc); a line con- 
nects the <V> in this plot. The halo stars with < Z < 
4 kpc (<Z> = 2.8 kpc) have <V> = -242+22 km s" 1 . 
This is compatible with them having (on average) zero ro- 
tation; it is those with 4<Z<8 kpc that have an excess of 
retrograde orbits. Fig. |4ja) and (b) include only stars with 
Z < 8 kpc and are to be compared with stars in the same 
Z-range shown in Fig. 1 (c) and (b) respectively in Majew- 
ski, Munn & Hawley (1996)(MMH) which are in the SA 57 
at the NOP. The vertical dotted line in Fig. |4]» and MMH 
(c) is the demarcation between prograde and retrograde or- 
bits. Both plots show a number of very retrograde orbits 
with negative W velocities. The pronounced correlation be- 
tween U and W in MMH (b) is present to a lesser extent 
in Fig. |4jb) 0. There is therefore qualitative agreement with 
MMH that there is structure in the phase-space distribution 
of these NGP stars. We test this quantitatively in the next 
Section. 

3.3 Is our NGP halo sample homogeneous? 

Table Ogives <U> and <V> and their corrected dispersions 
a u and o v for our 78 halo stars with Z< 8 kpc. Separate 
solutions are given for the RRL and BHB stars and also 
for the stars with positive and negative W space-velocities. 
These show that: 

• (a) The combined solution for <V> for both RRL and 
BHB stars and all W (—268+14 km s _1 ) is in good agree- 
ment with that found by Majewski (1992) (-275 km s" 1 ). 
The <V> of the BHB stars alone (-239+15 km s _1 ) is 
only mildly retrograde and compatible with zero halo rota- 
tion. The RRL sample alone, however, is strongly retrograde 

6 The different sign convention for U used by MMH causes nega- 
tive U to correlate with negative W in their case and positive W 
in our case. 
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Table 7. The RR stars of our sample. The U,V and W velocities (in km s 1 ) are heliocentric. The velocities Vr, V^> and V z (in km 
s _1 ) are in galactocentric cylindrical coordinates. For further details, see text. 



ID 


Z (kpc) 


R„„; (kpc) 

(J Lit. \ L J 


U 


V 


W 






v 2 


(1) 


(2) 


(3) 


('1) 


(5) 


(6) 


(7) 


(8) 


(9) 


GR-Com 


13.9 


17.3 


-035±095 


-470±111 




-001±095 


-239±111 




GH-Com 


8.9 


13.0 


+113±076 


-373±096 




+133±076 


-131±096 




IQ-Com 


10.4 


14.2 


-144±086 


-137±065 




-142±086 


+084±065 




NSV-5476 


7.4 


12.0 


+011±062 


-332±051 


-019±050 


+019±062 


-100±051 


+012±050 


V-Com 


3.4 


9.2 


-221±086 


-133±153 


-018±028 


-215±086 


+092±153 


+011±028 


CD-Corn 


9.4 


13.4 


-233±068 


-703±085 


-257±014 


-214±068 


-475±085 


+250±014 


AF-031 


9.5 


13.4 


-133±072 


-444±124 


-202±050 


-118±072 


-215±124 


+195±050 


TU-Com 


4.4 


9.7 


+194±095 


-368±168 


-074±016 


+204±095 


-134±168 


+067±016 


CK-Com 


6.6 


11.3 


+047±042 


-557±049 


-070±012 


+051±042 


-326±049 


+063±012 


AF-042 


14.1 


17.3 


-099±151 


-363±141 




-083±151 


-135±141 




CL-Com 


7.7 


11.9 


+076±083 


-207±050 


-025±016 


+085±083 


+026±050 


+018±016 


AT-CVn 


6.8 


11.4 


-089±086 


-266±121 


+089±052 


-081±086 


-032±121 


-096±052 


GY-Com 


11.4 


14.8 


+293±089 


-204±095 




+300±089 


+044±095 




GS-Com 


14.8 


17.6 


+111±186 


-335±293 




+126±187 


-095±293 




RR-CVn 


2.6 


8.8 


-053±031 


-390±040 


+014±011 


-047±031 


-157±040 


-021±011 


S-Com 


1.7 


8.3 


-055±034 


-205±033 


-066±004 


-046±034 


+027±033 


+059±004 


SV-CVn 


2.7 


8.8 


+090±018 


-326±028 


+081±028 


+095±018 


-098±028 


-088±028 


FV-Com 


(5.5 


10.7 


+180±060 


-154±061 


-171±050 


+190±060 


+076±061 


+164±050 


U-Com 


1.7 


8.2 


-264±020 


-314±019 


-038±003 


-255±020 


-083±019 


+031±003 


SW-CVn 


2.8 


8.8 


-157±036 


-200±036 


-012±021 


-146±036 


+038±036 


+005±021 


DV-Com 


7.3 


11.0 


+280±066 


-316±086 


-127±010 


+289±066 


-083±086 


+120±010 


AF-791 


(5.9 


11.0 


-362±061 


-036±054 


-230±012 


-334±061 


+226±054 


+223±012 


AW-Com 


10.5 


13.4 


+235±072 


-145±092 




+246±072 


+080±092 




TX-Com 


5.6 


10.0 


+236±052 


-357±048 


-088±049 


+240±052 


-134±048 


+081±049 


AP-CVn 


4.6 


9.5 


-181±038 


-247±044 


-009±049 


-173±038 


-007±044 


+002±049 


EM-Com 


9.5 


12.6 


+152±102 


-334±087 


-105±010 


+155±102 


-111±087 


+098±010 


AF-155 


7.8 


11.5 


+018±083 


-262±090 


-098±032 


+025±083 


-032±090 


+091±032 


TY-CVn 


4.0 


9.1 


-359±048 


-303±058 


+084±016 


-353±048 


-057±058 


-091±016 


IP-Corn 


7.2 


10.8 


+157±087 


-471±100 


-111±030 


+155±087 


-246±100 


+104±030 


SA57-19 


6.4 


10.2 


-110±049 


-364±052 


+059±050 


-104±049 


-130±052 


-066±050 


EO-Com 


(5.9 


10.6 


-154±078 


-459±064 


+089±010 


-153±078 


-222±064 


-096±010 


UV-Com 


10.1 


12.9 


-293±067 


-314±084 


+076±011 


-290±067 


-059±084 


-083±011 


TZ-CVn 


5.3 


9.7 


+165±075 


-361±074 


-131±050 


+166±075 


-139±074 


+124±050 


SA57-47 


5.9 


9.8 


+137±054 


-701±067 


-134±037 


+128±054 


-474±067 


+127±037 


SA57-60 


5.4 


9.1 


-020±058 


-473±068 


-090±050 


-017±058 


-241±068 


+083±050 


EW-Com 


8.0 


11.2 


-220±063 


-239±074 


+050±012 


-211±063 


+015±074 


-057±012 


IS-Com 


4.4 


8.9 


-047±032 


-440±044 


-132±010 


-045±032 


-207±044 


+125±010 


AF-882 


6.3 


10.3 


+043±047 


-286±050 


-071±050 


+044±047 


-061±050 


+064±050 



(-328±28 km s" 1 ). The dispersion a v of the 20 RRL stars 
with negative W (136±21 km s _1 ) is much greater than that 
(13±4 km s _1 ) of the 6 RRL that have positive W. Despite 
the rather small numbers, a variance ratio test shows this 
difference to be significant at better than the 99% level. The 
distributions of V for the BHB and RRL are shown in Fig. [5] 
The hatched histograms are for the stars that have negative 
W- velocities. Following Shapiro & Wilks (1965), separate 
tests of the V-distributions of the BHB and RRL stars with 
Z < 8 kpc show that neither have significant differences from 
normal distributions. 

• (b) There is an asymmetry in the ratio of RRL to BHB 
stars in the sense that the ratio is larger in the sample that 
has negative W velocities. This effect is significant at better 
than the 99% level and suggests that we are not dealing with 
a homeogeneous halo population. 

• (c) The <U> space-velocity is significantly negative 
(-72±21 km s _1 ) for both BHB and RRL stars that have 



positive W velocities. Those with negative W have the ex- 
pected zero <U>. 

There is therefore strong evidence that the NGP halo 
stars with Z< 8 kpc do not belong to a homogeneous popu- 
lation. It does not seem likely that observational selection in 
the samples could produce the differences in the kinematic 
properties between the BHB and RRL stars noted in (b) 
above. The RRL stars cover roughly the same area of the 
sky as the BHB stars (Fig. 1) and the six RRL with positive 
W are distributed widely in both R.A. and Declination. The 
lack of homogeneity must therefore be present over much of 
the survey area. We noted in Sec. 2.3 that stars with R.A. 
> 13 h :00 m showed predominantly negative radial velocities. 
There are 6 RRL stars with these R.A. and Z < 8 kpc. Their 
heliocentric space-velocities and dispersions are: 

< U > = +010 ±63 km s" 1 a u = 128 ± 37 km s' 1 
<V> = -417 ±74 km s' 1 a v = 152 ± 44 km s' 1 
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Table 8. Mean heliocentric space-velocities U, V & W (km s _1 ) for our program BHB & RRL stars as a function of height Z (kpc) 
above the galactic plane and for local samples of these stars. Local A is a sample of BHB stars (Kinman et al. 2000) and Local B and 
C are the HALOl and HAL02 samples of RRL stars (Martin & Morrison, 1998). Local D is the HAL02 sample after normalization to 
our distance scale. 



Range <Z> n„„ <U> a u <V> 

in Z 

(1) (2) (3) (4) (5) (6) 



(7) (8) 



<W> 
(9) 



(10) 



0-4 



4-8 



8 -12 



>12 



6.0 



14.1 



26 



52 



9.8 25 



14 



-52±23 110±15 
(-55±20) 

-12±20 131±13 
(-7±16) 



+1±41 

(+7±36) 

-26±40 
(-24±36) 



-242±22 103±14 26 
(-239±20) 

-282±18 lllill 52 
(-278±14) 



172±24 



328±34 119±17 21 
(-328±27) 



-349±24 
(-350±23) 



11 



±14 
(-1±20) 



0±33 
(-2 ±24) 



67± 9 



-37±13 85± 8 
(-37±11) 

-32±24 106±15 
(-28±22) 



96±22 



Local A 27 -1±25 129±18 -210±15 79±11 27 

Local B 81 -8±20 180±14 -197±12 111±9 81 

Local C 84 +1±21 193±15 -219±10 91±7 84 

Local D 84 +3±23 214±17 -240±11 100±8 84 



-14±20 101±14 

-8±10 93±7 

-5±10 96±7 

-8±11 108±8 



Table 9. Mean heliocentric space- velocities <U> & <V> and the total space-velocity <T> (km s x ) for our program BHB & RRL 
stars for different ranges of the space-velocity W. 



Star 


Range 




<U> 




<V> 




<T> 




Type 


in W 
















(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


RRL 


W< 


20 


+12±41 


166±26 


-321±36 


136± 21 


403±29 


125± 19 


BHB 


W< 


24 


+11±20 


78±11 


-216±24 


85± 12 


272±22 


104± 15 


All 


W< 


44 


+12±21 


124±13 


-264±22 


128± 14 


332±20 


131± 14 


RRL 


W> 


6 


-112±41 


135±39 


-351±31 


13± 4 


410±31 


69± 20 


BHB 


W> 


28 


-64±20 


108±14 


-258±19 


82± 11 


309±19 


96± 13 


All 


W> 


34 


-72±21 


112±14 


-275±17 


81± 10 


327±17 


99± 12 


RRL 


All W 


26 


-17±35 


165±23 


-328±28 


118± 16 


405±23 


113± 16 


BHB 


All W 


52 


-29±16 


102±10 


-239±15 


93± 9 


292±14 


101± 10 


All 


All W 


78 


-25±16 


125±10 


-268±14 


109± 9 


329±13 


117± 9 


RRL(P>0. d 6) 


All W 


8 


-8±66 


163±41 


_414±75 


189± 47 


490±50 


132± 33 


RRL(P<0. d 6) 


All W 


18 


-21±44 


171±28 


-290±21 


33± 6 


374±20 


82± 14 



< w > 



-085 ± 32 km s~ 



■ 059 ± 17 km s~ 



This small group thus shows a highly retrograde <V> 
and a negative < W >, but their corrected velocity disper- 
sions arc not unusual. 



3.4 The RRL sample 

The 26 RRL stars with Z< 8 kpc have a mean height above 
the plane of 5.3 kpc. Our analysis of the QSO proper mo- 
tions gave an rms systematic error in the proper motions 
of 0.3 mas yr _1 in each coordinate. This corresponds to a 
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Figure 4. Plots of V^el vs W/, e i and Uh e ; vs W^ e ; are shown in (a) and (b) respectively for the stars with Z < 8 kpc. The vertical 
dotted line in (a) divides stars with retrograde orbits (on the left) from those with prograde orbits. The lower figure (c) plots Vh e ; vs the 
height (Z) above the galactic plane for the whole sample. The locations of the mean V^el of Local A and Local D are shown by crosses. 
Filled triangles and open circles represent the RRL and BHB program stars respectively; the line connects the values of <V> for these 
stars. 
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Figure 5. The distributions of V for (a) the RRL stars and (b) 
the BHB stars in our NGP sample. Stars with hatched histograms 
are those with negative W velocities. 



~2 km s _1 error in V at a distance of 1 kpc. Consequently 
for our Z<8 kpc sample, a 2a systematic error in <V> from 
the proper motions should on average produce a ~20 km s _1 
systematic error in V. It does not seem likely therefore that 
systematic errors in the proper motions can wholly account 
for the retrograde motion (~100 km s _1 ) of the RRL sam- 
ple. The retrograde V could also be produced if our RRL 
distances are 45% too large. This corresponds to a 0.8 mag. 
error in the modulus of these stars. We used the same My 
for the RRL stars as for the reddest BHB starfl It does not 
seem likely, therefore, that our My can have such a large 
error. Neither does it seem likely that the mean apparent 
magnitudes of the RRL stars could have as much systematic 
error since they are based on multiple photoelectric obser- 
vations. 

The total space velocities (T) of our program stars were 
calculated from the space velocities U, V & W (Tables[S]and 
by: 



T = ^/U 2 + V 2 + W 2 



(8) 



The 15 BHB stars that have Z<8 kpc and (B - V)>0.10 have 



<V> = -227±35 and tr v 
retrograde motion. 



116±21 km s 1 so do not show a 
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Figure 6. Plots of the total space velocity (T) in km s -1 as a 
function of the distance from the galactic plane (Z) in kpc for (a) 
BHB stars and (b) RRL stars. The BHB stars at the NGP are 
shown by open circles and local BHB stars by filled squares. The 
RRL stars are shown by filled triangles; those with W-velocitics 
< —50 km s~ x are shown by encircled triangles. All stars above 
the line would be classified as belonging to the HAL02 sample of 
Martin & Morrison (1998). 



These total space velocities are plotted against the 
height above the galactic plane (Z) for the BHB stars and 
RRL stars in Figs. [6ja) and |6jb) respectively. The Local A 
sample of nearby BHB stars is also plotted as filled squares 
in Fig. E|a). The encircled triangles are the RRL stars that 
have W- velocities < —50 km s _1 . The line: 



T{kms~ 1 ) = 235 - 86 x Z 



(9) 



shows the boundary that is used by Martin & Morrison to 
separate their HAL02 from their DISK2 populations; all our 
program stars belong to their HAL02 population. Most of 
our stars have values of T that are well below the Galactic es- 
cape velocity of from 500 to 600 km s _1 (Smith et al. 2006). 
Possible exceptions are CK COM and SA 57-47 that have 
T = 563449 and 727466 km s _1 respectively. Large total 
space-motions are associated with (a) stars with 4<Z<8 kpc; 
(b) RRL stars rather than BHB stars; (c) RRL stars with 
large negative W-velocities and (d) RRL stars whose periods 
are greater than 0.60 days (Fig. [7] and Table [5J). Among the 
longer period RRL stars, the ones with lower amplitudes 
have the larger T (Fig. [7}. For periods greater than 0.60 
days, we would expect those of lower amplitude to belong 
to Oosterhoff Class I and the higher amplitude variables to 
belong to Oosterhoff Class II. The mean period < P a b > 
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Figure 7. The total space velocity (T) in km s — 1 for the RRL 
stars with Z<8 kpc plotted against the logarithm of the period P 
(days). The encircled triangles are of type c. A period of 0.60 days 
corresponds to log P = —0.221. The numbers by the symbols rep- 
resenting the longer period variables are their visual amplitudes. 



of the 20 type ab RRL that have Z < 8 kpc is 0. d 57. Our 
assumed RRL absolute magnitudes are linked to the RRL 
in the LMC whose < P ab > are 0. d 583 and 0. d 573 from the 
MACHO (Alcock et al. 1996) and OGLE (Soszyhski et al. 
2003) surveys respectively. The < P ab > of the RRL stars 
in Galactic globular clusters is 0. d 585 (Clement et al. 2001); 
65% of these variables are Oosterhoff type I 0. This sug- 
gests that the LMC contains predominantly Oosterhoff type 
I variables. This is confirmed by their Fourier types (Al- 
cock et al., 2004) and so our distances should be correct for 
Oosterhoff I variables which (judging from their < P ab >) 
probably comprise most of our sample. If, however, the Oost- 
erhoff II variables are as much as ~0.2 mag brighter than 
the Oosterhoff I variables (Lee and Carney 1999) then their 
distances and their total space velocities will have been un- 
derestimated. This may explain some of the dispersion in T 
for the RRL stars with P > 0.60 days. 



3.5 The Galactocentric space motions 

The Galactocentric space motions of our program stars in 
the cylindrical vectors Vr, and V z have been calculated 
according to the equations given by Kepley et al. (2006) 
which were chosen to be compatible with those used by 
Helmi et al. (1999). These velocities are given in cols. (7), 
(8) & (9) respectively of Tables [S] and [7J The Z-coordinate in 
this galactocentric system is in the opposite direction to that 
in our heliocentric system so that the galactocentric V z has 
the opposite sign to W in the heliocentric system. Galacto- 
centric V^, and heliocentric V are respectively negative and 
increasingly negative for retrograde orbits. The mean val- 
ues of these space- velocities <Vh>, <V^> & <V 2 > for the 



8 The RRL stars in most of the intrinsically faint, metal-poor 
dSph galaxies have < P ab > 0. d 60; the Sgr dSph on the other 
hand, has < P ab > = 0. d 574 (Cseresnjes 2003). 
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program stars with Z< 8 kpc and their corresponding dis- 
persions (after correction as described in Sec. 3.1) are given 
in Table [10] It is seen that the BHB stars show no galacto- 
centric rotation <V^> = — 005±015 km s _1 in agreement 
with the analysis of the radial velocities of 1170 BHB stars 
from the Sloan Digital Sky Survey (SDSS) by Sirko et al. 
(2004b). Table [lOl also gives the dispersions in cylindrical 
coordinates derived by Sirko et al. for three of their samples 
which we call SDSS(a), SDSS(b) and SDSS(c). These refer 
to their whole sample, stars with g < 18 and stars with g 
< 16 respectively. The dispersions for our BHB are similar 
to those found for the BHB samples of Sirko et al. (2004b) 
and have or ~ o^ whereas non-BHB samples of the local 
halo ellipsoid (Sirko et al. 2004b, Table 1) are like our RRL 
sample and have or > o>- There are only 6 RRL stars in 
our sample with positive V^, compared with 20 with negative 
velocities; the corresponding numbers for the BHB sample 
are 25 and 27. A Fisher's 2x2 test shows that this difference 
is significant at the 96%-level. A Shapiro- Wilks test of the 
V</>-distribution of the RRL stars shows no significant de- 
parture from normality. It therefore seems unlikely that the 
retrograde < V^> of these stars is caused by the inclusion of 
a few RRL stars with highly retrograde orbits but rather is a 
property of the whole sample. 

3.6 Stars with Z > 8 kpc 

In this paper we have chosen not to analyze the data for the 
39 stars with Z > 8 kpc. At 10 kpc, the random errors in 
the proper motions produce errors in the space motions that 
are comparable with the space motions themselves and the 
estimated 2a systematic error will be ~40 km s -1 . Currently 
we do not have radial velocities for all these fainter stars and 
those whose velocities have errors ^25 km s should be re- 
observed. 



4 SUMMARY AND CONCLUSIONS 

The kinematics of the BHB and RRL stars in the NGP show 
that this halo sample is not homogeneous. The BHB stars 
show essentially no galactic rotation in agreement with Sirko 
et al. (2004b) while the RRL stars show a definite retrograde 
rotation. Our whole sample of 117 stars includes stars out to 
16 kpc, but we restricted our discussion to the 52 BHB stars 
and 26 RRL stars with Z < 8 kpc whose space motions are 
the more reliable. This sample of 78 stars has a heliocentric 
<V> = — 268±14 which is in reasonable agreement with 
the retrograde rotation of —55 km s _1 found by Majewski 
(1992) and Majewski et al. (1996) for their SA-57 sample 
of subdwarfs at a similar <Z>. Both our BHB and RRL 
sample and their subdwarf sample show qualitatively similar 
dependencies of \Jhei on Whei and an excess of stars with 
negative Whei (Fig- HI- We also find that our BHB stars 
with Z< 8 kpc (<Z> = 4.7 kpc) have o r = 101±10, o> 
= 93±09 and o z = 81±08. These dispersions are similar to 
those of the BHB sample of Sirko et al. (2004b) with g < 18 
which has ay = 97±10, o> = 109±18 and o z = 101±6. These 
dispersions are more isotropic than local samples of halo 
stars that typically have o r >o^>o z . Qualitatively, this is 
what we might expect from a halo that is more flattened 
towards the plane and more spherical away from the plane 



and at greater galactocentric distances (Kinman et al. 1965; 
Wesselink 1987; Preston et al. 1991; Vivas & Zinn 2006). The 
local BHB sample (Local A in Table|5} has oy = 129±18, o^ 
— 80±11 and o z — 101±14 so that or > o<t, but (unusually) 
also has o z > o$. Many of the stars in this sample were 
selected as early-type high-velocity stars by Stetson (1991) 
and so Local A may have some kinematic bias. 

Our RRL stars that have positive W-velocities not only 
have a much smaller velocity dispersion (oy) than those with 
negative W-velocity, but the ratio of RRL to BHB stars is 
significantly lower in this group compared with the group 
that has negative W-velocities. This inhomogeneity there- 
fore involves a streaming motion that covers a significant 
part of our survey area and which is more pronounced for 
the stars with an RA > 13 h :00. A similar streaming occurs 
in the subdwarf sample of Majewski et al. (1996) in the 
same part of the sky. A model of the Sgr stream (Fig. 1 of 
Martmez-Delgardo et al. 2006) also predicts an infall in the 
Northern Galactic Cap but not quite in the direction of our 
survey and probably refers to more distant stars than most 
of those in our sample. The streaming that we have observed 
may be connected with the Virgo overdensity (Lupton et al. 
2005) but this overdensity is not yet sufficiently well defined 
for this to be more than a speculation. 

The simultaneous use of two halo tracers (BHB and 
RRL stars) has shown that the halo has two field star com- 
ponents that seem to parallel the old and young components 
of the halo globular clusters. This is not at all a new idea. 
Thus, Wilhelm et al. (1996) (who give references to previous 
suggestions for two-component halos) analysed the radial 
velocities of 525 BHB stars and found an "inner" flattened 
halo (Z<4 kpc) in prograde (+40 ±17 km s^ 1 ) rotation and 
an "outer" (Z>4 kpc) more spherical halo with a retrograde 
(-93±36 km s" 1 ) rotation . This relatively high retrograde 
motion for the "outer" BHB stars does not agree with either 
our results or those of Sirko et al. (2004b). A Z>4 kpc sam- 
ple will be at high galactic latitudes where the stars lie far 
from the apices of galactic rotation and their radial velocities 
are relatively insensitive to the galactic rotation vector. This 
insensitivity of the radial velocities and also possible inho- 
mogeneities in the halo may be responsible for this discrep- 
ancy. On the other hand, the analysis of the space motions 
of halo stars by Chiba & Beers (2000) found no evidence for 
a retrograde "high halo" ; possibly because their sample did 
not extend far enough from the solar neighbourhood. 

The two halo components that fit our data are similar 
to those discussed by Borkova and Marsakov (2003) from 
their analysis of the space-motions of nearby RRL stars. 
One component is associated with a strong blue horizontal 
branch, has no galactic rotation and has a relatively isotropic 
velocity dispersion. The other component is characterized 
by having RRL stars whose < P a b > shows that they are 
mostly of Oosterhoff type I. This second component has, on 
average, a retrograde rotation and a less isotropic velocity 
dispersion than the first component. Our second component 
also has streaming in its W space-motion which suggests 
that it consists of tidal debris. This does not mean that the 
first component contains no RRL stars or the second has 
no BHB stars but that the BHB stars predominate in the 
first component and the RRL stars in the second. We would 
expect that the first component, having a blue HB, would 
be associated with Oosterhoff type II RRL stars. We must 
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Table 10. Mean galactocentric space-velocities Vr, & V z (km s 1 ) and their corrected dispersions an, cr^ & c z for our program 
BHB & RRL stars and for comparison the SDSS BHB stars of Sirko ct al. (2004b) (a) whole sample, (b) with g < 18, and (c) with g < 
16. 



Star n <V_r> a R <V > cr <V Z > a z 

Type 

(1) (2) (3) (4) (5) (6) (7) (8) 



RRL 


26 


-11±35 


162±23 


-095±29 


122± 17 


+043±17 


079± 11 


BHB 


52 


-20±16 


101±10 


-005±15 


93± 9 


+005±12 


081± 08 


All 


78 


-17±15 


124±10 


-035±14 


111± 9 


+017±10 


081± 09 



SDSS(a) 1170 
SDSS(b) 773 
SDSS(c) 227 



94±08 
97±10 
109±20 



098± 18 
109± 18 
122± 24 



098± 05 
101± 06 
094± 12 



be careful about this assertion, however, because there is 
considerable scatter in the HB-type vs. < P a t > plot; e.g. 
M 62 has < P ab > = 0. d 548 (Oosterhoff type I) and yet has 
a blue HB (Clement et al. 2001; Contreras et al. 2005). 

It would be very desirable to be able to compare all 
the properties ([Fe/H], space densities, HB type, Oosterhoff 
type etc) of the halo stars in the North Polar Cap with 
those in the solar neighbourhhod. This is beyond the scope 
of this paper and would require considerable attention to 
selection effects. As an example, all RRL surveys miss vari- 
ables with amplitudes below a certain limit. Our photome- 
try of the BHB stars showed that a number of these were 
low-amplitude RRL stars so, in this respect, our polar sam- 
ple is unusually complete and contains low-amplitude long- 
period type ab and type c variables that were missed in 
earlier surveys. Our results do not necessarily conflict with 
those of earlier surveys that have sampled larger volumes of 
space. The volume of space that we have sampled is rela- 
tively small. At the mean distance of 5 kpc, the width of the 
field is no more than ~1 kpc. This may account for the pro- 
nounced lack of uniformity of the halo that we have found at 
the NGP and makes it essential that similar surveys should 
be made at other key locations such as the South Galac- 
tic Pole and the Anticenter and if possible that they should 
include more than one halo tracer. 
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Table 11. The BHB stars of our sample. The coordinates are for J2000. The parallaxes TIhbp, ^hbc an d TIhbk are defined in Sect. 
2.2.3; the error ajj refers to the adopted parallax TIhbp- Proper motions and their errors (columns 12-15) are in masyr^ 1 . Radial 
velocities and their errors (columns 16 and 17) are in km s~ 1 . 



n_sr.BC 


^HVK 


Ma 




M<5 




RV 


o-RV 


(10) 


(ii) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


0.280 


U.2oO 


—15.8 


1.9 


1 o 
— t.i) 


1.1 


A Ct 

4b 


50 


0.067 


0.067 


1.0 


A 1 

4. 1 


—3.8 


a a 
4.0 


— to 


a n 
4U 


0.332 


0.282 


—9.2 


3.5 


—20.0 


1.5 


28 


10 


0.100 


0.090 


2.4 


3.4 


—13.7 


6.2 


—92 


10 


0.064 


0.086 


—3.8 


2.2 


—4.0 


1.5 


128 


40 


0.231 


0.228 


1.4 


1.6 


—2.0 


1.8 


—24 


10 


0.081 


0.074 


-1.9 


2.1 


-5.9 


1.0 


—40 


10 


0.123 


0.121 


—3.6 


1.9 


—7.8 


1.9 


—32 


a n 
4U 


0.150 


0.144 


U.4 


2.8 


—3.9 


i i 
1. ( 


—62 


10 


0.178 


U.152 


1 Q 

— ( .5 


O A 

z.4 


—y.y 


1.0 


A O 

42 


10 


n i 11 
U.l / / 


U. 104 


A Q 

— 4.o 


1.0 


n a 

— y.4 


2.0 


on 


OU 


0.120 


0.104 


—6.8 


i i 
1. / 


i n 

— / .y 


3.3 


1Q 

— to 


10 


n nQi 
U.Uo ( 


U.UOz 


—o.y 


1.0 


— l.U 


1.0 


Oo 


a n 
4U 


0.134 


0.132 


—9.2 


1.6 


—9.0 


1.9 


95 


10 


0.296 


0.281 


—4.9 


1.2 


—2.0 


1.3 


1 1 
1 7 


10 


0.091 


0.092 


—3.9 


1.6 


—2.0 


3.4 


64 


10 


0.301 


0.248 


-7.8 


2.1 


—9.9 


1.8 


22 


10 


0.065 


0.066 


—2.0 


3.2 


A O 

— 4.y 


2.2 


2 t 


4U 


U.l i 2 


n i rzi 
U.15 t 


—o.y 


A O 

4.z 


—6.9 


O 1 

2. ( 


etc. 
— 00 


10 


0.153 


0.137 


—1.8 


1.3 


—13.0 


2.9 


11 


10 


0.083 


0.083 


a n 
— 4.U 


1.8 


—3.9 


1.0 


—103 


10 


0.333 


0.278 


—0.5 


I i 
1. / 


A O 

—4.2 


1.3 


—85 


10 


U. L2o 


U. IZo 


— o.U 


1 Q 
l.O 


— O.U 


O.O 


—10 


a n 
4U 


0.133 


0.112 


—4.9 


2.5 


—6.0 


3.1 


—120 


a n 
4U 


0.173 


0.163 


—4.0 


2.6 


—5.8 


1.5 


—134 


10 


0.151 


0.135 


—5.5 


1.4 


—4.8 


2.1 


181 


10 


0.060 


0.058 


1.0 


1.6 


-4.0 


1.5 


—41 


40 


0.122 


0.116 


0.3 


o o 
z.o 


—10.9 


2.0 


—50 


1U 


0.102 


0.097 


—1.8 


2.5 


—6.0 


n i 


07 
O I 


10 


0.110 


0.103 


—0.6 


3.1 


—8.2 


a n 
4.U 


—159 


10 


0.230 


0.211 


2.1 


0.9 


—5.9 


2.0 


—178 


10 


0.293 


0.262 


—17.8 


2.0 


—16.6 


2.5 


9 


10 


U.o ( t 


n qi a 


— 10. ( 


O Q 
2.6 


— iy.o 


Q Q 
O.O 


on 

— zy 


1U 


0.263 


0.247 


—16.9 


1.5 


—9.1 


2.2 


t 


10 


n i ic\ 
U. 1 ( U 


U. loO 


i o 

— / .y 


O Q 
2.0 


—0.2 


1 A 
1.4 


"I i 
1 1 


An 
4U 


0.277 


0.245 


-17.9 


1.8 


-20.9 


1.8 


-70 


10 


0.203 


0.195 


-12.0 


1.9 


3.2 


1.9 


-18 


10 


0.406 


0.340 


1.0 


1.9 


-9.0 


1.6 


2 


10 


0.614 


0.598 


-25.9 


1.9 


-35.8 


1.7 


4 


10 


0.127 


0.115 


-5.8 


1.7 


-8.6 


1.2 


13 


10 


0.348 


0.324 


-5.9 


2.5 


-3.8 


4.1 


145 


10 


0.088 


0.087 


-1.0 


1.9 


-4.0 


0.9 


106 


10 


0.109 


0.104 


-4.0 
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47 


10 
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ID 



RA 



DEC 



V 



B-V E(B-V) 



K 



^HBP 



(1) 




(2) 






(1) 




(4) 


(5) 


(6) 


(7) 




(8) 


l 


16549-51 


12 


01 


02 


.1 


34 


37 


40 


13.67 


0.078 


0.016 


13.30 


A 


0.271 


0. 


AF-003 


12 


03 


20 


.1 


32 


01 


38 


16.48 


0.230 


0.019 


15.35 


C 


0.066 


0. 


AF-006 


12 


04 


13 


1 


29 


41 


50 


13.64 


0.007 


0.016 


13.54 


A 


0.309 


0. 


AF-011 


12 


06 


19 


,5 


31 


51 


37 


15.85 


0.101 


0.019 


14.72 


U 


0.098 


0. 


AF-013 


12 


07 


11 


,2 


28 


58 


52 


16.78 


0.115 


0.019 


17.18 


U 


0.063 


0. 


AF-727 


12 


08 


12 


.0 


36 


02 


08 


13.93 


0.137 


0.019 


13.52 


A 


0.229 


0. 


AF-022 


12 


10 


08 


,2 


30 


53 


00 


16.37 


0.082 


0.019 


15.96 


U 


0.078 


0. 


AF-729 


12 


10 


26 


.4 


34 


32 


54 


15.24 


0.146 


0.013 


14.70 


A 


0.123 


0. 


AF-029 


12 


12 


59 


,8 


33 


43 


13 


14.85 


0.131 


0.013 


14.22 


A 


0.148 


0. 


AF-030 


12 


13 


01 


,9 


29 


09 


05 


14.85 


0.044 


0.023 


14.48 


A 


0.169 


0. 


16022-26 


12 


11 


29 


.2 


27 


23 


52 


14.65 


0.097 


0.023 


14.18 


A 


0.172 


0. 


AF-038 


12 


15 


01 


,6 


32 


02 


21 


15.63 


0.049 


0.016 


15.23 


B 


0.114 


0. 


AF-039 


12 


15 


01 


.8 


30 


06 


18 


16.65 


-0.006 


0.019 


15.66 


D 


0.081 


0. 


AF-041 


12 


15 


35 


.4 


32 


35 


38 


15.02 


0.168 


0.013 


14.26 


A 


0.134 


0. 


AF-045 


12 


16 


19 


,9 


28 


56 


06 


13.45 


0.129 


0.026 


12.86 


A 


0.292 


0. 


AF-046 


12 


16 


53 


,7 


30 


04 


09 


15.98 


0.130 


0.023 


15.69 


D 


0.090 


0. 


AF-048 


12 


17 


20 


9 


32 


11 


27 


13.92 


-0.008 


0.013 


13.85 


A 


0.278 


0. 


AF-050 


12 


17 


34 


,8 


29 


19 


24 


16.57 


0.222 


0.026 


15.89 
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0.065 


0. 
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12 


18 


05 
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17 


25 
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0.081 


0.019 


14.30 
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37 
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0.026 
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22 
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36 
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0.019 


16.19 
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0. 
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12 
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02 


.6 


27 


27 


16 


13.65 


0.009 


0.019 


13.50 


A 


0.310 


0. 
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12 
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24 
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36 


31 


09 
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0.013 


14.84 


A 


0.126 


0. 
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12 


23 


26 
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36 


41 


12 


15.46 
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0.013 


15.06 
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0.126 


0. 
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12 


26 


10 
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30 


07 


13 


14.63 


0.108 


0.019 


14.12 
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0.170 


0. 


AF-070 


12 


26 


37 
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28 


49 


21 


15.19 


0.047 


0.023 
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0.144 
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12 


28 


56 
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32 
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29 
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0.016 
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0.059 


0. 
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30 


33 
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30 
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15.36 


0.115 


0.016 


14.89 
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0.120 


0. 
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11 
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30 


29 


29 


15.72 
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0.019 


15.03 


B 


0.101 


0. 
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51 
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32 


31 


07 
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0.013 


15.10 


B 


0.108 


0. 
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12 


31 


10 


.6 


30 


51 


26 


14.05 


0.092 


0.016 


13.51 


A 


0.224 


0. 
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12 


31 
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57 


01 
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0.062 


0.019 


13.32 


A 


0.281 


0. 
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12 


36 


21 
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27 


16 


34 


13.40 


0.006 


0.019 


13.28 
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0.350 


0. 
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37 


53 


13.74 


0.090 


0.013 


13.35 


A 


0.257 


0. 
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12 
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.7 


33 


52 


27 


15.22 


-0.014 


0.016 


15.09 


C 


0.157 


0. 
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12 


41 


17 
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29 


09 


15 


13.86 


0.035 


0.013 


13.67 


A 


0.262 


0. 
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12 


41 


38 
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33 


20 


55 


14.33 


0.087 


0.016 


14.08 


A 


0.198 


0. 
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12 


42 


07. 


9 


30 


55 


41 


13.21 


0.006 


0.016 


13.09 


A 


0.378 


0. 


AF-113 


12 


42 


12 


,4 


31 


56 


48 


11.78 


0.139 


0.016 


11.22 


A 


0.611 


0. 


AF-115 


12 


42 


39 


3 


33 


23 


52 


15.42 


0.068 


0.016 


15.02 


B 


0.123 


0. 


16466-15 


12 


43 


41 


,8 


32 


21 


28 


13.14 


0.089 


0.013 


12.71 


A 


0.339 


0. 


AF-130 


12 


11 


18 


,6 


29 


26 


20 


15.92 


0.180 


0.013 


15.25 


B 


0.088 


0. 


AF-131 


12 


11 


18 


.6 


30 


55 


22 


15.58 


0.119 


0.016 


15.01 


B 


0.108 


0. 


AF-134 


12 


11 


12 


,2 


32 


29 


55 


12.76 


0.170 


0.013 


12.04 


A 


0.378 


0. 



on 
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Table 11. continued 



ID RA DEC V B-V E(B-V) K n HBP cr n n HBC n HVK n a an a n 6 crfig RV crRV 



(1) 




(2 


) 




(3) 




(4) 


(5) 


(6) 


(7) 




(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


16031-44 


12 


46 


15.6 


27 


10 


47 


13.80 


0.145 


0.013 


13.23 


A 


0.238 


0.010 


0.239 


0.234 


1.4 


1.2 


-13.0 


3.2 


11 


10 


15622-48 


12 


46 


43.7 


27 


17 


17 


14.47 


0.153 


0.013 


13.97 


A 


0.174 


0.007 


0.174 


0.174 


-3.1 


1.4 


-10.2 


2.6 


-69 


10 


15622-07 


12 


46 


47.3 


27 


27 


44 


12.86 


0.121 


0.013 


12.44 


A 


0.374 


0.016 


0.378 


0.369 


-5.1 


2.3 


-12.9 


3.3 


-37 


10 


AF-138 


12 


47 


37.2 


33 


07 


28 


14.72 


-0.009 


0.013 


14.65 


A 


0.193 


0.017 


0.209 


0.171 


-8.9 


3.0 


-10.5 


1.6 


101 


10 


15622-09 


12 


48 


58.8 


26 


48 


04 


13.08 


0.122 


0.010 


12.48 


A 


0.336 


0.014 


0.339 


0.324 


-1.7 


1.4 


-14.0 


2.5 


-121 


10 


AF-797 


12 


50 


23.2 


36 


20 


53 


15.42 


0.074 


0.016 


15.07 


A 


0.122 


0.006 


0.126 


0.116 


-8.0 


2.2 


-6.0 


0.3 


85 


40 


AF-804 


12 


51 


52.7 


36 


40 


42 


14.07 


0.098 


0.016 


13.64 


A 


0.221 


0.010 


0.225 


0.212 


-13.0 


2.3 


-5.9 


1.0 


94 


10 


SA57-001 


12 


52 


18.7 


28 


19 


58 


14.43 


0.155 


0.010 


13.76 


A 


0.176 


0.007 


0.176 


0.173 


-1.9 


2.8 


-18.4 


1.8 


-124 


10 


SA57-006 


12 


53 


57.3 


29 


38 


08 


13.99 


0.083 


0.010 


13.58 


A 


0.229 


0.011 


0.235 


0.218 


-14.6 


2.7 


-5.6 


2.2 


112 


10 


SA57-007 


12 


54 


03.7 


28 


51 


47 


12.18 


-0.010 


0.013 


12.11 


A 


0.626 


0.058 


0.678 


0.553 


-13.9 


1.6 


-28.5 


2.2 


-92 


50 


SA57-009 


12 


54 


18.0 


28 


43 


42 


15.96 


0.204 


0.013 


15.14 


B 


0.085 


0.003 


0.085 


0.085 


-9.6 


2.7 


-9.1 


1.7 


47 


10 


SA57-017 


12 


56 


23.8 


27 


28 


39 


13.41 


0.084 


0.016 


13.00 


A 


0.304 


0.015 


0.313 


0.289 


-25.4 


4.7 


-10.6 


3.5 


-4 


10 


SA57-021 


12 


57 


34.6 


31 


48 


45 


16.62 


-0.066 


0.016 


16.43 


U 


0.097 


0.011 


0.107 


0.069 


-6.0 


2.2 


5.8 


3.3 


-70 


40 


SA57-029 


12 


59 


13.6 


30 


36 


19 


16.31 


0.085 


0.010 


15.72 


D 


0.079 


0.004 


0.081 


0.073 


-2.9 


1.6 


-4.0 


1.8 


-182 


40 


SA57-032 


12 


59 


43.3 


31 


20 


30 


15.13 


0.155 


0.013 


14.63 


A 


0.128 


0.005 


0.128 


0.128 


-7.8 


1.5 


-8.0 


2.5 


-14 


10 


SA57-036 


13 


00 


41.2 


28 


52 


37 


16.23 


0.062 


0.010 


15.37 


B 


0.084 


0.005 


0.087 


0.075 


-3.9 


1.9 


-5.9 


2.7 


-5 


10 


SA57-040 


13 


03 


13.6 


28 


40 


11 


15.14 


0.080 


0.010 


14.67 


A 


0.135 


0.007 


0.139 


0.127 


-2.9 


1.5 


-2.5 


1.2 


30 


10 


SA57-041 


13 


03 


20.6 


28 


43 


28 


16.33 


0.092 


0.010 


15.62 


C 


0.077 


0.004 


0.079 


0.072 


-4.9 


1.9 


-3.5 


1.8 


127 


10 


SA57-045 


13 


04 


22.7 


31 


30 


25 


15.31 


0.108 


0.010 


15.12 


B 


0.122 


0.005 


0.123 


0.125 


1.0 


1.4 


-7.9 


1.0 


-135 


10 


AF-825 


13 


.04 


56.3 


35 


03 


31 


14.32 


0.051 


0.013 


14.08 


A 


0.207 


0.013 


0.216 


0.196 


-5.9 


0.8 


-5.1 


1.4 


-139 


50 


SA57-046 


13 


05 


08.7 


27 


48 


01 


16.71 


0.148 


0.013 


16.76 


C 


0.062 


0.002 


0.062 


0.072 


-4.9 


1.7 


-4.4 


2.2 


151 


40 


SA57-055 


13 


07 


19.4 


30 


25 


56 


16.03 


0.140 


0.013 


15.58 


C 


0.086 


0.003 


0.086 


0.085 


-1.0 


3.1 


-5.9 


0.7 


-224 


10 


AF-841 


13 


08 


30.3 


32 


43 


40 


14.06 


0.091 


0.013 


13.62 


A 


0.222 


0.010 


0.227 


0.212 


-7.8 


1.2 


-10.1 


1.7 


-25 


10 


AF-848 


13 


10 


48.5 


32 


28 


19 


11.75 


0.001 


0.010 


11.64 


A 


0.732 


0.061 


0.786 


0.661 


-9.8 


3.0 


-18.0 


2.8 


-116 


10 


SA57-066 


13 


11 


29.4 


31 


47 


06 


14.68 


0.116 


0.013 


14.08 


A 


0.163 


0.007 


0.165 


0.156 


-3.0 


3.4 


-3.9 


2.1 


-243 


10 


AF-854 


13 


12 


01.6 


37 


26 


55 


15.11 


0.058 


0.013 


14.77 


A 


0.142 


0.008 


0.148 


0.133 


-10.9 


1.8 


-0.0 


1.4 


-10 


50 


SA57-080 


13 


13 


10.2 


30 


45 


00 


14.16 


0.077 


0.013 


13.90 


A 


0.215 


0.011 


0.222 


0.210 


-8.0 


1.2 


-4.8 


2.5 


-46 


10 


AF-866 


13 


13 


39.3 


34 


27 


18 


13.93 


0.059 


0.010 


13.64 


A 


0.243 


0.014 


0.252 


0.230 


-6.0 


1.8 


-11.0 


2.8 


-41 


50 


SA57-087 


13 


14 


29.3 


31 


42 


15 


15.97 


0.093 


0.010 


15.45 


C 


0.091 


0.004 


0.093 


0.086 


-4.5 


2.6 


-3.8 


2.6 


-112 


10 


AF-900 


13 


20 


25.1 


33 


57 


56 


13.71 


0.013 


0.013 


13.55 


A 


0.293 


0.023 


0.313 


0.265 


-9.6 


1.5 


-10.8 


0.6 


92 


50 


SA57-111 


13 


20 


37.2 


30 


12 


14 


16.47 


0.094 


0.016 


15.66 


U 


0.073 


0.003 


0.075 


0.068 


-5.9 


1.3 


-4.0 


1.5 


-138 


10 


AF-909 


13 


22 


12.0 


31 


21 


31 


15.26 


0.023 


0.013 


15.09 


A 


0.141 


0.010 


0.149 


0.130 


-1.9 


1.3 


-5.9 


1.2 


-144 


10 


AF-914 


13 


23 


21.0 


37 


47 


05 


12.27 


-0.003 


0.010 


12.16 


A 


0.581 


0.050 


0.626 


0.519 


-36.2 


3.4 


-15.0 


3.4 


95 


50 


AF-916 


13 


2d 


15.4 


33 


22 


13 


14.00 


0.181 


0.010 


13.69 


A 


0.210 


0.008 


0.211 


0.222 


-0.2 


1.3 


-6.0 


1.7 


-132 


50 


AF-918 


13 


25 


02.6 


29 


46 


18 


14.08 


0.028 


0.013 


13.87 


A 


0.240 


0.017 


0.254 


0.220 


-7.9 


2.3 


-5.9 


1.6 


-24 


50 



Table 12. The RRL stars of our sample. <K> magnitudes are corrected for phase effect as described in the text. The parallaxes labelled 
ITy, life and TIhb were estimated as described in Sect. 2.2.2. The error an refers to the adopted parallax Tly. RV and crRV are in km 
s" 1 . 



ID RA DEC logP [Fe/H] <V> <K> E(B-V) Ii v a n n K n HB fi a a^ a fig a/x s RV crRV 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) 



GR-Com 


12 


05 


18.6 


27 


56 


58 


-0.282 





16.33 


15.14 


B 


0.021 


0.071 


0.004 


0.075 


0.072 


-3.8 


1.3 


-6.0 


1.7 








GH-Com 


12 


05 


24.5 


27 


38 


22 


-0.352 





15.38 


14.42 


A 


0.020 


0.110 


0.007 


0.114 


0.112 


-1.7 


1.5 


-8.9 


2.3 








IQ-Com 


12 


06 


03.6 


27 


59 


18 


-0.294 





15.70 
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